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^ ■ Abstract. The C function of a massive scalar field near a cosmic string is computed 

I and then employed to find the vacuum fluctuation of the field. The vacuum expec- 

tation value of the energy-momentum tensor is also computed using a point-splitting 
approach. The obtained results could be useful also for the case of self-interacting 
' scalar fields and for the finite-temperature Rindler space theory. 
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I. Introduction 

In the last few years many authors have considered quantum fields on spacetimes contain- 
ing conical singularities. These kind of manifolds are relevant for studying fields in presence 
of real conical singularities, namely existing in the Lorentzian section of the manifold, as in 
the case of idealized cosmic strings and of 2-|-l-dimensional gravity. In other cases, conical 
singularities appear in the finite temperature theory, when adopting the Euclidean path inte- 
gral with periodic imaginary time formalism. It is the case of spacetimes containing horizons, 
such as the Rindler and the Schwarzschild spaces, which are much studied in the off-shell 
approach to the one-loop quantum corrections to the Bekenstein-Hawking entropy. It is 
worth noticing that, in these latter cases, the conical singularity could be avoided adopting 
a canonical approach. 

As a consequence of this interest in conical spaces, many techniques have been developed 
to compute relevant quantities in presence of conical singularities (see, among the others, 
i, I, I, I, i g 0, I, I, 0, 0). With few exceptions @, |, H, [l, HI, 0, HI, most authors 



considered massless fields only. Indeed, the introduction of the mass complicates considerably 
the problem, so that it becomes very difficult to obtain manageable forms for the physical 
quantities. At this regard, it has to be noted that using the heat-kernel approach M, M 
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it seems that all the problems vanish: in fact, for a scalar field of mass m the integrated 
heat-kernel is simply related to the massless one as 



—tm? T^m=0 



(1) 



and, at least for a simple cone ^ x Ca, the integrated massless heat-kernel is well known 
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where 27r — a is the deficit angle of the cone, V{Ca) is the volume of the cone and S£)-2 is 
the volume of the transverse dimensions. On the contrary of the massless case, the Mellin 
transform of Eq. (|l|) exists in the massive case, and so it is easy to compute the massive 
global C, function and from it quantities like the effective action, which have a simple closed 



form. However, it is well known |11] that this procedure does not yield the correct dependence 
on the deficit angle of the cone in the massless limit and, moreover, these global quantities 
are finite, apart for the volume divergences, while the local quantities show a non-integrable 
singularity near the tip of the cone.0 Therefore, it is at least dubious that this simple result 
for the massive case is correct. Besides, in order to compute important quantities such as the 
vacuum fluctuation of the field and the stress tensor, local rather than global quantities are 
needed. 

The aim of this work is to give a more manageable tool to compute local quantities for 
a massive scalar field around a cosmic string, and we will accomplish this by computing 
the local C function, which will be given as an expansion in powers of mr, where r is the 
distance form the string core. Therefore, the obtained expression is useful near the string, 
where mr <C 1.0 However, since the local quantities diverge at the string core, this is just 
the interesting region if, for instance, one wants to consider the back-reaction of the energy- 
momentum tensor on the background metric 0, 18]. The expression we obtain is simple 
enough to allow us to compute the renormalized vacuum expectation value {(fP'{x)) and the 
the energy-momentum tensor, potentially up to arbitrary order in (mr)'^. Notice that the 
most common tool used to compute local quantities is the Green's function, which in the 
massive case is given in a complicated integral form |12, 13, ll^, 0] or as a sum of generalized 



hypergeometric functions |16| 



The capability of treating the massive case is useful also when dealing with self-interacting 
theories, such as Xcj)'^. In fact, consider an Euclidean path- integral approach to the theory of 
a self-interacting scalar field: the generating functional of the theory is the given by 



Vcj)expi-S[<P]), 



where V{4>) is a potential. In the usual 0^-theory, it has the form V{(l)) = ^m^cl)'^ + ^A(/>^, 
A being the coupling constant. In the one-loop perturbative approach, the above action is 
expanded near the minimum (p of the action up to second order terms 

1 



^In the cosmic string case, this problem could be overcome renormalizing the action of the string as in 
^Note that r = h/mc is the Compton wavelength of the particle. When the result are translated for the 
Rindler case, the approximation is valid close to the event horizon. 
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where A is the Laplace-Beltrami operator and = V"{(p). In the (p'^-case, we have that 
= m? + iA</)^. Then, the standard step is to assume (p is a. constant configuration, so 
that standard passages lead to effective action 

W = - \nZ = S[4)] + ^ lndet[-A + M^] + 0{n^), 

where, in general, S'[</>] is proportional to the volume of the spacetime. The determinant in the 
second term is ultraviolet divergent and it has to be defined by means of some regularization 
procedure. Convenient ways to regularize the above determinant are, for example, the heat- 
kernel plus Schwinger's proper-time regularization or the (^-function regularization. 

From the above outline, it is clear that the one-loop self-interacting case is closely related 
to the massive case, and, in particular, the regularization of the determinant is the same, 
provided that m? is replaced by M^. Of course, the renormalization of the effective action 
to remove the ultraviolet divergences is different. In this paper we shall not discuss the self- 
interacting case and its renormalization explicitly, but we observe that in the cosmic string 
case, in which {(lP'{x)) oc r^^, it could be more reliable a 'Hartree-like' approximation as that 
employed in ||l^ and in which in the action A(/>^ is replaced by \{(jP')(jP' . 

As we said above, we are mainly concerned in the comic-string background. We remind 
that the space-time around an infinitely long, static and straight cosmic string [p!s| ] has the 
topology = E? X Ca, where Ca is the simple two-dimensional cone with deficit angle 
2ir — a, and the metric is 

The polar angle deficit 27r — q is related to the mass per unit length of the string fj, by 
27r — a = SirGfi. Throughout this paper we shall assume positive deficit angle, so that 
u = In /a > 1. For GUT strings — 1 ~ 10~^ In the rest of the paper we adopt an 

Euclidean approach and so we perform a Wick rotation t = it: the form of the metric is the 
same as above with the replacement —dt'^ — > ^dr'^. Moreover, for technical reasons we will 
also consider the obvious generalization to the D-dimensional case, with topology R^~^ x Ca- 

We have assumed zero thickness for the string: this is clearly an idealization, since an 
actual cosmic string should have a small but finite radius, of the order of the Compton wave- 
length of the Higgs boson involved in the phase transition which gives rise to the cosmic string 
fl^ . The internal structure of the string has non-negligible effects even at large distances, 
as it has been shown by Allen, Kay, and Ottewill However, in the same work it has 

been shown that for a minimally coupled scalar this dependence on the internal structure 
is absent, and one can safely use the idealized string. In order to avoid the complications 



discussed in |2C], which are related to the non-uniqueness of the self-adjoint extension of the 
Laplace-Beltrami operator in the idealized conical space-time, we shall consider the mini- 
mally coupled case only. Actually, in Sec. IV we will compute the stress tensor and, for sake 
of completeness, we will consider arbitrary coupling. Therefore, we have to remember that 
for ^ 7^ others effects could be present in a realistic cosmic string. 

Although the results we will obtain are easily translated to the case of fields at finite 
temperature in the Rindler space Q or near the horizon of a Schwarzschild black hole, we 

*The Euclidean metric of the cosmic string and of the finite-temperature Rindler space are in fact the same. 
In this latter case, 6 corresponds to the periodic Rindler time and the period a is the inverse temperature. 
The absence of the conical singularity fixes the Unruh-Hawking temperature T — 1/2tv in suitable units ^]. 
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prefer to consider the cosmic string background only. This because the Euchdean Rindler 
space or Schwarzschild black hole show a conical singularity only when the temperature is 
different from the Unruh-Hawking one, namely when one considers the 'off-shell' theories, but 
it has been shown that the off-shell quantum states of a field are affected by several pathologies 



on the event horizon |22] and that these theories do not have a consistent thermodynamics 



The rest of this work is organized as follows. In Sec. II we compute the C function of a 
massive scalar field in the cosmic string background. Then we use this result to compute, in 
the region mr <C 1, the vacuum fluctuation of the field in Sec. Ill and the expectation value 
of the stress tensor in Sec. IV. Section V contains the conclusions. In the rest of the the 
paper we use units for which h = c = G = 1. 

II. Computation on the ( function 

We consider a quasi-free real scalar field in the background given by the (Wick-rotated) 
metric . The action of the theory is the given by 

= J d^xV5^^ 

where A is known as the small fluctuation operator and in our case reads 

A = -Ao + m'^ + (,R. (3) 

Here ^ is a parameter which fixes the coupling of the field to the gravity by means of the scalar 
curvature R. We shall consider the minimally coupled case, ^ = 0. Then, by means standard 
passages one sees that the generating functional of the theory can be expressed in terms of the 
functional determinant of the small fluctuations operator, which can be conveniently defined 
by means of the C function regularization |24|: 

InZa = -^lndetn~^A = ^C'iO\ Afi'"^) 

where ^(•sl^) is the global C function related to the operator A, the prime indicates the 
derivative with respect to s and fi is an arbitrary parameter with the dimensions of a mass 
needed for dimensional reasons and not to be confused with the mass per unit length of the 
string. The global C function can be formally written as the integral over the manifold of a 
local C function ({s\A){x): 

Cis\A) = [ as\A)ix)^d^x. 
Jms 

Actually, when the manifold is non-compact only the local C function has a precise mathemat- 
ical meaning, since the integration requires the introduction of cutoffs or smearing functions 
to avoid divergences. 

In the massless case, the ( function on the cone has been explicitly computed in [|ll] ]. In 
order to compute it in the massive case, the starting point is the well known relation between 
the local ( function and the heat kernel given by the Mellin transform: 



1 
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The heat kernel of the massive field is related to the massless one by means of the following 
obvious decomposition: 

Now we take into account the following property of the Mellin transform of the product of 
two functions [25| 



roo 1 ra—ioo 

/ f-^f{t)g{t)dt = — / F{z)G{s - z)dz, 



where F and G are the Mellin transforms of / and g respectively and o" is a real number in 
the common strip of convergence of the two Mellin transforms. In this way we can write the 
massive C, function of in terms of the massless one: 

1 ra—ioo 

ris; x) = —— / r(z)r=°(z; x)rr?^-^'T{s - z)dz, (4) 
where we have used the Mellin transform 



/■oo 

/ t^-^e"*"^ dt = m-2T(s), 
Jo 



which converges for Res > 0. 

Since we are interested in the Q function of a massive scalar field in the (Euclidean) space 
]^D~2 ^ (j^^ consider the corresponding massless C, function. It has been shown in Q and 
pA[| that it is the sum of two parts which converge in separate strips of the complex s-plane 
and which are summed only after their analytic continuation: 

r=0(z;x)=C<(5;x) + C>(5;x), 



where 



r 



2s-D 



. , . ■ r(s-^)r(^-s) D-\ D 

(47r)°— r(s) 2^{s-i^) 2 2 

^2s-D Yis - D D 

C>(s;x) = ^ ' Go^{s-^), <Res<- + z., 

a(47r) — r(s) 2 2 

and the function Ga{s) is defined as 



n=l 



and can be analytically continued in the whole complex plane showing simple poles in s = 1, 
with residue a/47r, and s = Un + k + l, k = 0,1,2, ... , (if a 7^ 27r) with obvious residue. Since 
C< and C> do not have a common strip of convergence, we must split Eq. (^ in two parts: 
setting Res > D/2 we have 

C^is;x) = r-^m-^^ r dzr{z--^)G^{z--^ns-z){mrr 

2mr[s) (47r) — 0Fa Jf<Rez<Res 
2s-D ( 00 (_i\n 
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where we have performed the integral shifting the integration path to the right and picking 
up the residues of the poles of the integrand. The same procedure can be applied to C<^(*; x), 
but in this case we have to consider the poles in z = s + n and z = 2 + n: 



^<("''^^ - 2vrzr(.)(4^)^^«yi-l<Re.<Ref 2r(z + 1 - Z)/2) '^^"^'^ 

i g (-ir („,,)2nr(. + n - -s-n) 



(4^)^0Far(s) IS) 2T{s + n + l-D/2) 



n=0 



n! ' ' 2r(n + l) 



Now we analytically continue each term and sum to get the final expression of the local C, 
function of a massive scalar field: 

(47r) — ar(s) 
(mrr-2s g (-ir ^^^^2„r(n + l/2)T{s - n - D/2) 



n n\ 2r(n + l) 



n=0 

^1^0 ^' r(2i.„ + A; + l) J 

(5) 

where the function Ia{s) is defined as 

r(i-«) 



Us) = 



vr 



G„(s) + 



2r(«) 



The function is analytic in the whole complex plane but in s = 1, where it has a simple 
pole. Near the pole we have 



where 7 is the Euler's constant. In expression (^) we have also reintroduced the arbitrary 
mass fi using the formal relation = M^^C(s|^)) aiid which has been omitted in 

the derivation for simplicity. It must be noted that the above C function can be obtained in 
a more direct way, albeit much longer, performing the integrations in the Mellin transform 
of the spectral representation of the massive heat kernel, and then rearranging the sums 
and the Euler gamma functions in the generalized hyper geometric functions obtained by 
the integrations. We prefer the used method since it can be more easily applied to the 
computation of others quantities, as we will see in the next sections. 

Although expression ^ looks awful, it is very simple when mr <C 1: considering the 
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physically interesting case D = 4 and terms up to0 (mr)^ 



47rar(s 
87rQr(s) 



/«(s - 1) - (mr)2/^(s) + ^(mr)4/^(s + 1) + O ((mr)") 



+ 



r(. 



{mry 



r(s-3) 



+ 0{{mrY 



In the first term of the first row we recognize the massless C, function on B? xCq, |11|, while the 
first term in the second row becomes the C function of a massive scalar field in the Minkowski 



space-time when a = 27r [26|. The others terms are clearly corrections due to the presence 
of the conical singularity. When the singularity is absent, namely when a = 27r, the first row 
vanishes, since I2-k{s) = 0, while in the second and third rows only the first term survives, 
since the others cancel two by two. 

At the same order in mr, the effective lagrangian density is given by 



C{x) 



2ds 



L_|/^(_i)_(^^)2j^(o) 



(mr) 
2~ 

(mr) 
2A 



— If rn\ \ni> m 
7 + ln^ + +ln — 

V \ I ) V /X 



2Ga(2) + 27 + 21n 



mr 
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+ 2(mr 



^/ir(2i/ + i) 



where Gain), n > 2, are readily computed being in the region of convergence of the series 
defining Ga- For example 



GU2) 
G„(3) 



1 

2^ 
1 



27 + ^(z^-^) 



where ipix) = r'(x)/r(x). In the limit a ^ 2tt the above effective lagrangian reduces to the 
usual Coleman- Weinberg potential. It is important to note that the computation of higher 
orders in the above expansion in powers of (mr)^ does not involve particular complications. 



III. Vacuum fluctuations 



Now we use the above C function to compute the (renormalized) value of the vacuum 
fluctuation of the field p^ : 

((/)^(x)) = limC(s;x). 

Since even in the Minkowski space the vacuum expectation value diverges, we renormalize 
the vacuum expectation value on the cone subtracting the Minkowski value: in this way we 

^Note that we are interested in the values of the ^ function at s — 0, 1, and we consider positive deficit 
angles only, so that u > 1. 
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have a remarkable cancellation of the poles, yielding a finite result: 



1 



[mr 



TflT 1 



mr 



2G„(2) + 27-l + 21n — 



(7) 



The computation of higher corrections is straightforward. One can verify that Eq. (0) 
correctly vanish in the limit a — > 27r, and that up to order (mr)^ it is identical to the result 



obtained by Moreira Jnr [16|. In particular, we notice the additional logarithmic divergences 
of the vacuum fluctuations at the conical singularity due to the massive corrections. 



IV. Energy-momentum tensor 

Another important vacuum average is that of the energy-momentum tensor. While in 
the massless case it is quite easy to be computed, since its form is fixed from symmetry 
arguments, in the massive case the computation is much more difficult and, to our knowledge, 
the massive corrections near the string have never been shown explicitly. Only in |13| and 
1 15] the explicit form of the energy-momentum tensor has been given for mr ^> 1, where they 
found the expected exponential damping factor exp(— 2mr). 

We use a point-splitting approach, in which the vacuum expectation value of the energy 
momentum tensor is given by the coincidence limit of a non-local differential operator applied 
to the propagator of the field ||2l| : 



{T^uix)) =i lim Df^i^{x,x')G:F{x,x') 



where 



D^,{x,x') = (1 - 20V^V,, + -2CV^V, + (2e - 1/2)5, 



m 



(8) 



and the prime indicates that the derivative has to be taken with respect to x' . Since the 
propagator can be obtained from the off-diagonal C function as 

Gjrix, x') = i lim C,{s; x, x'), 



we can compute the energy momentum tensor from the C function as (see also |24, p^j) 

{T^u{x)) = - lim lim D^^{x, x')C{s; x, x'). 



The partial derivatives of the off-diagonal C function which appear in the above expression 
can be computed from the spectral representation in D dimensions (x = (t, z), k = |k|) 



Cd{s;x,x') 



2(47r) 



D-2 

2 



ar( 



D-2^ 



dk k 
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2n , 



dX X\X'^ + k"^ + m^] 
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and then taking the coincidence hmit. In this way one can easily show that 

ded0/CD{s;x,x')\x=x' = -9|Cd(s; x, a;')U=x' 
dz,d^,'CD{s;x,x')\^=^' = -d^^CD{s;x,x')\a:=x' = 27rCD+2{s;x) 

dTdr'CD{s;X,x')\x=a:' = -drCois] X , x')\x=a:' = '^■^Cd+2{s; x) 

dr'CD{s;x,x')\x=x' = -drCD{s;x). 
As far as dQ(^j:){s; x, x')\x=x' is concerned, in the m = case it is easy to see that we have 

^6»Cd is;x,x)\x=x'= Ha{s — ^f2.^^ 

{A7r)^~ ^/7^ aT{s) 

where the function Ha{s) is defined and studied in the Appendix A. The massive case can 
then be treated using the off-diagonal version of Eq. (^) with the partial coincidence limit 
r = r', z = z' and t = t': 

dlUs;e,9')\e=e' = / r(z)a|cr=0(z; 0, 0')|e=e'r(s - z)m2^-2^d^ 

2mT[s) jRc z>3 47rVvrQ 

2s-4 2s f oo (_-,^n 

= - + " - 3/2)fl„(. + n - 1) 

hh r(2.„ + fc + l) I 

where, as usual, we have performed the integration shifting the integration contour to the 
right and picking up the residues at z = s + n and z = Un + k + 2 to get an expansion in 
powers of mr similar to Eq. (^). 

As far as the second derivatives with respect to r and r' are concerned, using the following 



identity, which can be proved using some recursion formulas for the Bessel functions |25 



one can see that 



r zr 



1 „ „ . , . 1 



drdr'C,D{s;x,x')\x=x' = —drrdrCois; x) + -^dgCD{s;x,x')\x=x' +Xd{s 



2r 



X 



1 „ . . . 1 



dr'CDis;X,x')\x=x' = - — drCDis;x) KdgCD{s]X,x')\x=x' - XDis;x) 

zr r 



where the function xd{s] x) is defined in Appendix B. If D = 4 and m = we simply have 

xS=4°(«;^) = 4vr(.-2)G=0(s;x), (9) 
while the massive case is more complicate and is studied in the Appendix. 
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Now we have all the pieces needed to compute the massive correction to the energy- 
momentum tensor. For the actual calculation it is convenient to follow |29] and |30|. We 
define the renormalized stress tensor as 

which, considering only the first massive correction, turns out to be 



(Trrix))^ 

(Tttix))^ 



1 



47rar^ 

-1 
47rar^ 

-1 

(Tuix))^, 



2H^{0) - (6^ - l)/„(0) + (mrf H^{1) + ^ 



2i?,(0) - - {2i - 1)7,(0) + {mr f H^{1) + ^ 



U-l) + 2(4e - l)/„(0) + (mr)2l,(0) 



(V(^)) 



R 



The tensor can be written in the more familiar form 

-1 
liio^rM 

+lb{mrf{v - 1)(12^ - 1 - z^)diag(0, 1, -1, 0) 
-lb{mrf - l) diag(l, 0, 0, 1) 



(v^ - l) diag(l, 1, -3, 1) + 10(6^ - 1) _ 1^ diag(2, -1,3, 2) 



(10) 



In the limit m ^ the result is in agreement with that obtained by others authors ||3|, ^] 
The components {Trr{x))^ and {T0q{x))^ satisfy the equation 



{rT\)=T%, 

ar 

which follows from the conservation law V^T'^j^ = 0. It is interesting to note that in the 
corrections to the stress tensor of order (mr)^ are not present the logarithmic divergences 
which are instead present in the vacuum fluctuations at the same order (see Eq. (0)). The 
logarithmic terms appear in higher order corrections, which can also be computed with the 
method we have developed. Actually only the logarithms in the term of order (mr)^ give rise 
to divergences, since at higher orders they are multiplied by positive powers of r. Indeed, the 
correction of order (mr)^ to {Tgg{x))^ is 



1 



47rar2 



^ ^ 2 ^ + (2^ - l)^(r) + 4e - 1 



+ {mrf+^''u[2C + (4^ - 



r(-z. - l)r(z. + 1/2) 
V^T{2u + 1) 



from which one can also compute the correction to {T„{x))^ by means of the conservation 
law. The correction to {Tu{x))^ = {Tzz{x))^ is 



1 



^[2i?(r)-(4e-l)(A(r) + l)-i 



+(mr) 



2+2^2 4e - ly ^ / \ i/ • 

^y^TT{2u + 1) J 



10 



In the above equations we have set 

TflT 

A{r) = 2G<,(2) + 27 + 21n— . 

„ , , V — 1 f mr 1 \ In z/ 

^^^^ = ^r^+^-2j+— ■ 

An interesting point to be discussed is the dependence on the parameter ^ which fixes 
the couphng of the scalar field with the gravity. In the introduction we said that we would 



consider the minimally coupled case only, ^ = 0. This because it has been shown in |2C] that 
the idealized conical space is a good model of the space time of a cosmic string only for the 
minimally coupled case. For nonminimally coupled fields quantities like {(fP'{x)) will depend 
on the details of the metric in the core of the string even very far from the string. However, 
one could be interested, e.g., in finite-temperature fields in the Rindler space, where there is 
a true conical singularity and not just an idealization of a non-singular metric. In these cases 
it is interesting to consider also nonminimally coupled fields. 

From the mathematical point of view, it is not clear which is the meaning of the field 
equation [—A + rin? + £,B\4' = when the curvature R has Dirac's delta singularities, and 
our choice of = allowed us to avoid the problem. A possible way to define the problem 
is to smooth out the singularity, as done in |5|, 10 1: as a result, one finds that when the 



regularization of the singularity is removed to recover the conical space also the dependence 



of the Green's function on the parameter ^ vanishes |3C]. Then we can argue that also the C, 
function is independent of ^, since the modes used to construct the C, function are essentially 
the same as those for the Green's function. 

The dependence on the parameter ^ comes back into play when considering the stress 
tensor. It is worth while noticing that even when the manifold if fiat everywhere (in our case 
there is a Dirac's delta singularity in the curvature at the core of the string) the parameter ^ 
remains in the theory as a rehc of the fact that Tjj^y is obtained by varying the metric g^u in 
the field Lagrangian |21]. One can then see that, if ii = 0, the global conserved quantities as 



total energy should not depend on the value of ^. This is because the contributions to those 
quantities due to ^ are discarded into boundary surface integrals which generally vanish. 
However, this is not the case dealing with the cosmic string because such integrals diverge at 
the origin. Notice that similar problems appear working in subregions of the Minkowski space 



in presence of boundary conditions |21|. In the point-splitting approach, the stress tensor is 
obtained applying the ^-dependent operator D^iy{x,x') (see Eq. (^)) to the Green's function 
or to the C function. Therefore, if the above argument that the Green's function (or the C 
function) is independent of holds, it is not contradictory to compute the energy-momentum 
tensor applying D^^^x, x') with ^ 7^ to the the Green's function computed setting ^ = 0, as 
we have done above and as done by most authors. 



V. Conclusions 



In this paper we have studied the C function of a massive scalar field in a cosmic string 
background, and we have obtained an expression, Eq. (P), which is useful in the region 
near the core of the string, mr <^ 1. By means of this expression we have computed the 
massive corrections to the vacuum fluctuations, Eq. (0), and to the energy-momentum tensor, 
Eq. (p^, up to order (mr)^, going beyond the known results. Higher corrections are also 
computable. 
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Possible extensions of this paper are the inclusion of a magnetic flux carried by the string, 
which gives Aharonov-Bohm effects, and the case of spin-1/2 fields. Also the limit mr ^ 1 
is worth studying, for example rewriting the C function (|5|) in terms of the hyper geometric 
function iF2[a; b, c; z] and then using its asymptotic behaviour for large z. However, the vac- 
uum fluctuations and the energy-momentum tensor in this limit have been already obtained 



by others authors with different methods [13, 15]. 



As a possible application of the results of this paper we point out the study of the back- 
reaction of the energy momentum tensor of a massive scalar field on the background metric 



of the cosmic string, as done by Hiscock |17| and Guimaraes |15]. The computation seems 
feasible since the first massive correction to the energy momentum tensor does not involve 
logarithmic terms. We hope to cover this topic in a future paper. 
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Appendix A 

In this Appendix we study the function Ha{s), which is defined as the analytic continu- 
ation of following series: 

It can be studied and analytically continued proceeding exactly as for the function Ga{s) 
studied in the Appendix of |11]. Then one sees that the series converges for Res > 2 and 
that the analytic continuation for [n/2] < Res < 2 (here [n/2] represents the integer part of 
n/2) is given by 

[n/2]+l 00 

c,(s)i.3-2-2,<.^(2s + 2j - 3) + 5] ylfn{yk, s), 

j=0 k=l 

where is the Riemann ( function and the function /„(z^, s) is generally unknown, but 
vanishes for s = 1/2,0,-1/2,-1,.... The coefficients Cj(s) vanish for s = —n/2, (n = 
— 1,0, 1,2,...) for all j > {n + l)/2, and the first ones have been given in |11]. 

The function Ha{s) has then a simple pole in s = 2 and near this pole we have 

i/.(.) = ^^ + ^(7-ln.)+0((.-2) 

Moreover, it has simple poles at s = 1 + Vn + k, k = 1,2, . . ., due to the gamma function 
in the numerator of the terms of the series, with obvious residue. Finally, it is possible to 
compute the value of the function Ha{s) for some useful value of s: 
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Appendix B 



In this Appendix we study the function xois', x), which is defined by the analytic contin- 
uation of the following spectral representation 

Xd{s;x)= '^yy,,, rdkk^-' f rdXX'[X' + k' + mr'Jl{Xr), (11) 
ar{{D - 2)/2) Jo Jo 

which is the same as C,Dis;x) with dAA — > dAA^. Note that the massless case is trivial, 
since in that case the function x is simply related to the C function, see Eq. (^. In the 
massive case, we can proceed in analogy to the massive function, namely employing Eq. 
considering D = A and using Eq. (|9|) 

1 ra— too 

x'BUs-.x) = TT-T^/ nz)x'Bzl{z-,x)m''^-^'T{s-z)dz 



2TiiT{s) Ja 
27rir{s) 



too 
a— too 



r{z){z - s)CB=6iz;x)m''-''T{s - z)dz. 



Then we go on as usual splitting Q'j^^q{z;x) as C = C< + C> find performing the integrals 
shifting the integration contour to the left and picking up the residues of the poles. The final 
result is an expansion in powers of mr: 

2s-6 r oo (_^\n 
^ ' Kn=0 

,^(-1)''/ , ^6+2n-2s^(n + l/2)^(s-n-3) 

^0 nl 2^r{n + 1) 

^ ' tito 2^T{2vn + k + l) j 

(12) 
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